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SUMMARY/ABSTRACT 


The report presents a definition of a VOR/DME, airborne and ground 
systems simulation model. This description was drafted in response to a need 
in the creation of an advanced concepts simulation in which flight station 
design for the 1980 era can be postulated and examined. The simulation model 
described herein provides a reasonable representation of VOR/DME station in 
the continental United States including area coverage by type and noise 
errors. The detail in which the model has been cast provides the interested 
researcher with a moderate fidelity level simulator tool for conducting 
research and evaluation of navigator algorithms. Assumptions made within the 
development are listed and place certain responsibilities (data bases, 
communication with other simulation modules, uniform round earth, etc.) upon 
the researcher. 


INTRODUCTION 

In order to better facilitate future flight station design in light of 
rapidly advancing technology, a joint effort involving two National 
Aeronautics and Space Administration centers is underway in conjunction with a 
major aircraft company. In this effort, a preliminary design of an advanced 
concepts simulator was formatted and is being constructed at each partici- 
pating organization. While individual research focuses will be maintained at 
each organization, a high degree of commonality has been a goal in the design 
so that maximum cooperation and less costly exchanges can more effectively 
occur during the research phases. 

The software and hardware development efforts for such a simulation 
facility are massive and thus the burden has been shared among the three 
organizations. The aircraft and ground systems software has been broken down 
into modules in order to maintain flexibility and growth. Specific module 
requirements and definition have been assigned to organizations. The 
requirements and definition of the VOR/DME systems module is the subject of 
this document. The requirements and definitions, including some limiting 
assumption, which are listed, are aimed at a moderate level of fidelity. 


LIST OF VARIABLES AND UNITS 


Input 



PVORi , 1=1 .2 

Logical , 

A/C power status, VOR 

PDMEi , i =1 ,2 

Logi cal , 

A/C power status, DME 

Lat (A/C) 

(deg) 

A/C Latitude 

Lon (A/C) 

(deg) 

A/C Longitude 


AlM-JQSO&F 



ALT (A/C) 

(ft) 

A/C altitude mean sea level 

Vg (A/C) 

(kts) 

A/C ground speed 

VORFi, i = 1 , 2 

Mhz 

Selected frequency 

DATA BASE all for VORFi, ■ 

i = 1 or 2 


VLATi 

(deg) 

station latitude 

FLONi 

(deg) 

station longitude 

ALTi 

(ft) 

station altitude, 


X 

mean sea level 

Identi 

(Text) 

call letters, morse code 

Typei 

0=T, 1 =L , 2=H 

station classification 

MVARi 

(deg) 

magnetic variation 

flagl i 

Logical 

station on/off status, VOR 

flag2i 

Internal 

Logi cal 

station on/off status, DME 

TBi 

(deg) 

True bearing (station to A/C) 

GRi 

(nm) 

ground range (great circle) 

SRi 

(nm) 

slant range 

C0i 

(deg) 

cone of confusion angle 

n v 

(deg) 

VOR noise term (total) 

N r 

(deg) 

VOR noise term (receiver) 

N t 

(deg) 

VOR noise term (transmitter) 

N c 

(deg) 

VOR noise term (course 
roughness) 

°r 

(deg) 

VOR standard deviation 
( recei ver) 


(deg) 

VOR standard deviation 
(transmitter) 

t 

( sec) 

time 

X 

(sec) 

time correlation parameters 
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°c 

(deg) 

VOR standard deviation (course 
roughness) 

w c 

rad/sec/kts 

VOR frequency (course 
roughness) 

w v 

dimension- 

less 

white noise source 

N p 

(nm) 

DME noise term (total) 

n b 

(nm) 

DME noise term (transmitter 
and receiver) 

a D 

(nm) 

DME standard deviation 

(transmitter and receiver) 

n dc 

(nm) 

DME noise term (course 
roughness) 

W DC 

(rad/sec/kts) 

DME frequency (course 
roughness) 

°DC 

(nm) 

DME standard deviation (course 
roughness) 

Output 



VORFi 

(Mhz) (To 
data base) 

selected frequency 

BOi 

(deg) 

mag. bearing (station to A/C) 

DMEOi 

(deg) 

slant range 

VFLOi 

logical 

validity status, VOR 

DFLOi 

logical 

validity status, DME 

R o 

ft 

Mean Earth Radius 
(20,887,749.4 ft) 


VOR/DME System Modeling Description 

This description describes a simulation model for station and receiver 
portions of a VOR/DME system. 

The advanced concepts simulator contains the provision for simultaneously 
tuning two VOR/DME stations. The information received from these ground 
stations is used by other modules within the simulation to manage navigation 
functions. The VOR and DME portions are contained in the same module but can 
be failed individually. This approach is consistent with the assumption that 
the VOR and DME occur in pairs only. The VOR/DME models are activated by the 
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selection of VHF frequencies on the integrated comm/nav management panel. 

This automatically serves as the selection process for the DME channel. The 
frequency ranges will generally be between 108 to 118 MHz. The output of the 
station will be a bearing and slant range from the station location of the 
inquiring aircraft. The update rate is normally 1/15 of a second. The 
station bearing will contain the magnetic variation associated with the 
station Jobation. Distances for the DME portion will be slant range in 
n miles (6076.1 ft/nm). 

The model assumes the existence of a data base containing relevant 
parameters Given selected frequency, the data base will provide station 
!' de '. } on ^l u ' de > a " d altitudes, as well as type, identification and an 
on/off failure flag. Station and aircraft locations are combined to determine 

irln ,t 3 h h ! range ;- The range ’ a1titude and station type are compared to 
IC/\0 standards for a valid region check, including the cone of confusion. No 

consideration of ground features (such as mountains, hills, buildings etc. 
that could occult the signals) is included in the modeling process. 

Noise and mean bias errors are added to the signal bearing and slant 
range Parameters governing the magnitudes will be suggested but should be 
alterable by any experimenter. The bias errors (both VOR and DME) will 

attempt to represent both station and receiver terms plus a colored course 
roughness term. 

There are two types of failures that can be introduced in this model- a 
power failure can disable either the VOR or DME portion of the model and is 
assumed to be passed on to the appropriate other modules, a station failure of 
® V0R ^ DME P rocess imbedded in the data base will cause a nonvalid 

J n -f S1mi ] ar ™ anner as an outside-proper-region coverage failure. The 
latter failure will be passed on in the validity variable. 


Summary of Assumptions and Conditions 
o VOR and DME stations occur together 

o Valid regions of coverage given by United States standard (station 

type) and no transient effect concerning in/out of coveraqe are 
modeled 3 

o Noise and mean bias are gaussian 

o Magnetic variation included 

o Provisions included for a/c power and/or station failure 
o Module suuported from appropriate data base. 

o Ground features and/or higher order error terms are not included. 
Figure 1 is a basic block diagram of the module. 
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Preliminary Design Layout 


A preliminary version of the ordering of events in the model is proposed 
in table 1. This process could be repeated for the second frequency or both 
frequencies handled simultaneously. The data base is assumed to exist and 
contains as a minimum the information as supplied in the check case setups of 
this paper. 


Logic Details 

The names assigned to variables in this document should be changed by the 
coders or systems integrators to properly correspond to the global variables 
structures. Hence, the names used herein are to be considered temporary for 
purposes of illustrations. The output logic variables will be validity flags, 
VFLOi , i = 1 ,2 and DMEFOi , i =1,2. The logic for hoth VOR and DME is split 
into two phases: (1) an A/C power status, and (2) station status or coverage. 

These are separated in that a data hase call can be avoided if A/C power to 
both the VOR and DME is in a failed state. Station status and regional 
coverage are treated deeper within the modules and must involve a data base 
cal 1 . 


Let VFLOi = status of valid tests for VOR including coverage, cone of 
confusion, or station failure i = 1,2; and 

DFLOi = status of valid tests for DME including coverage, cone of 
confusion or station failure i = 1 ,2. 

Then IF{ VFLOi or PVORi) then IDENT = no signal (Blank) and BOi - o + n v , 

I F ( DFLOi or PDMEi) then DMEOi = o + n D 

where PVORi = power status of VORi , determined via input logic i =1,2, and 
PDMEi = power status of DMEi , determined via input logic, i =1,2. Note it is 
assumed that power status logic to VOR/DME is routed elsewhere as required for 
caution and warning. 


Mathematical Equations 
For each VORFi i = 1 ,2. r 

GRi (ground range) = 60x COS" 1 SIN(VLATi (station) )SIN(Lat(A/C) ) 

+COS(VLATi)COS(Lat(A/C)COS(Lon(A/C)-VLONi)J 

TBi (true bearing) = COS -1 fsiN( Lat ( A/C) )C0S(GRi/60)SIN(VLATi ) " 

[ SIN (GRi/60)COS(VLATi ) 

IF SIN ( VLONi - Lon) < 0 then TBi = 360 - TBi. (If LONi and Lon(A/C) are 
given in negative values for west such as United States then this logic must 
be reversed) 

SR-j (slant range) = Abs [ (ALT(A/C)-ALTi (station) ) 2 

+ 4 SIN 2 (Ro+ALTi (station) )(Ro+ALT(A/C))] 1/2 

BOi (Rearing output) = TBi + n y + MVARi 
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The error term n y includes receiver, transmitter and a course roughness 
component such that n y = N p + N t + N c 

N p is a random constant, selected from a normal distribution, mean = zero and 
o r = 0.2 (deg) N p = N [0,o p ]. 

N t is likewise a random constant term such that N t = N [0,o t ] where 
Oj. = 0.15 (deg). 


The course roughness component is treated as colored noise with correla- 
tion time inversely proportional to A/C ground speed (V g ) 

E [N c (t)N c (t + t)] = a c 2 e " WcV gl T l 

let a c = 0.2 (deg) and w c = 4 x 10“ 3 rad/sec/kts. 

Note that the course roughness error is generated from Gaussian white noise 
W v (t) as shown below. 



and the autocorrelation function of the driving white noise is 
E [W v (t)W y (t +x)] =(2 o 2 c w c V g ) 6 (t ) 
where % = time correlation parameter. Correspondi ngly 


DMEOi = SRi + N d 


where 


N D = N B + N DC 

Ng = N [0, a Q ]for a D = 0.5 nm or Ng = ( . 03SRi ) 

whichever is greater, and N DC = is similar to N c of the V0R process with 
equivalent definitions for W D , and a DC . Specifically 
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W D (t) 



where the DME range error is modeled as colored noise, i.e., 

E pl DC (t)N DC (t + x)J = cr 0C e J-l- 

with t D = 400 sec. and o D q = 0.1 nm. 

As shown in the above diagram, the DME correlated error is generated from 
white noise with the autocorrelation given by E W_(t)W^(t + t) J= 6(x). 

For implementation guide see Appendix B. L u J 


Coverage Limits . 

The following is an excerpt from reference 1 concerning the VOR/DME 
coverage volumes. These volumes shall be used to define valid regions for 
this model. The station type or class designator will be obtained from the 
data base. 

Standard Service Volumes (SSV) .- Ground stations are classified according 
to their intended use. These stations are available for use within their 
service volume. Outside the service volume, reliable service may not be 
available. For standard use, the airspace boundaries are called standard 
service volumes. They are defined, in the table below, for the three station 
classes. 


SSV CLASS 
DESIGNATOR 

T (Terminal) 


L (Low Altitude) 


ALTITUDE AND RANGE BOUNDARIES 


From 1000 feet (305 m) AGL up to and including 

12.000 feet (3,658 m) AGL at radial distances out to 
25 n. mi. (46 km). See Figures 4 and 5. 

From 1000 feet (305 m) AGL up to and including 

18.000 feet (5,486 m) AGL at radial distances out to 
40 n. mi. (74 km). See Figures 3 and 6. 


H (High Altitude) From 1000 feet (305 m) AGL up to and including 

14,500 feet (4,420 m) AGL at radial distances out 
to 40 n. mi. (74 km). From 14,500 feet (4,420 m) 
AGL up to and including 60,000 feet (18,288 m) at 
radial distances out to 100 n. mi. (185 km). From 

18.000 feet (5,486 m) AGL up to and including 

45.000 feet (13,716 m) at radial distances out to 
130 n. mi. (241 km). See Figures 2 and 6. 


These SSV's are graphically shown in Figures 2 through 6. 
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Within 25 n. mi. (46 km), the bottom of the T service is defined by the 
curve in Figure 5. Within 40 n. mi. (74 km), bottoms of the L and H 
service volumes are defined by the curve in Figure 6. (Note metric measure- 
ments are given for convenience and are approximations.) The distance 
parameter to be compared against the defined boundaries is: 


qi = [R 0 + ALTi (station) + ALT(A/C)] SIN 


If qi is outside the defined boundary, then VFLOi and DFLOi = false. 

The following excerpt from reference 1 defines the cone of confusion to 
be included within this module. 

Vertical Angle Coverage Limitations .- Within the operational service 
volume of each station, azimuth signal information permitting satisfactory 
performance of airborne components is normally provided from the radio horizon 
up to an elevation angle of approximately 60° for VOR components and 
approximately 40° for TACAN components. At higher elevation angles, the 
azimuth signal information may not be usable. Distance information provided 
by DME will permit satisfactory performance of airborne components from the 
radio horizon up to an elevation angle of 60°. Thus given CO = 60° then 
VFLOi and DFLOi = false, where COS(CQi) = qi / SR i . In addition, a ground 
station failure can be inserted in the data base through the FLAGli and FLAG2i 
variables such that 

If FLAGli = false then VFLOi = false. 

If FLAG2i = false then DFLOi = false. 


Test Conditions 

The following cases are designed to exercise the V0R/DME module by 
specifying the input conditions and examining the expected outputs. Data base 
values are also supplied for appropriate conditions. Results are contained in 


Objective of test 

Power failure logic, bearing quadrants and noise sources 
Remaining quadrants of bearing and various ranges 
Valid coverage regions for station types 
Station failure 

In some of the above cases, noise sources are eliminated in order to unmask 
possible sources of error. In other cases, the noise sources themselves 
provide a proper indication and should be examined for the bounds. Tables 2 
and 3 contain the necessary data base and aircraft information for checks. 


Appendix A. 
Case # 

1 . a,b 

2 . 

3. a-c 

4. a,b 


GRi 

2Ro 


(see Fig. 7) 
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Concluding Remarks 


The simulation model described herein provides a reasonable representa- 
tion of VOR/DME stations in the continental United States. The detail in 
which the model has been cast provides the interested researcher with a 
moderate fidelity level simulator tool for conducting research and evaluation 
of navigator algorithms. Assumptions made within the development are 
consistent with other portions of the Advanced Concepts Simulation and place 
certain responsibilities (data bases, communication with other modules, etc.) 
upon the researcher. 


REFERENCES 

1. U. S. National Aviation Standard for the VOR/DME/TACAN SYSTEMS. 
September 2, 1982, Dept, of Transportation Federal Aviation 
Administration, 9848.1. 
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TABLE 1.- PRELIMINARY DESIGN FLOW 

Proposed step by step process: 

(a) Check master power/proceed/exit 

(b) Compare frequency to previous/proceed/skip (optional) 

(c) Issue call to data base subroutine 

(d) ^flag f MVAR data ^ ^ ^ l0 " 9 ’ aU ’ type ’ a " d indent > 

(e) Check flag/proceed/issue non-valid 

(f) Calculate bearing and range 

(g) Determine validity; given range, type, ALT/Exit 

(h) Add bias and noise (function of type) 

(i) Add MVAR (station location, data base) 

(j) Output, valid discrete, ident code, slant range, bearing 
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TABLE 2.- PSUEDO BATA BASE 


Flat Rock VOR/DME 

IDENT = FAK, MVAR = 6.5W 

113.3 (ch80) frequency 

LAT = 37° 31 min 30 sec (37.525°) 

LON = 77° 49 min 30 sec (77.825°) 

ALT = 400 ft 

TYPE = H, FLAG1 , FLAG2 

CAPE CHARLES VOR/DME 

IDENT = CCV, MVAR = 8°W 

112.2 (ch50) frequency 

LAT = 37°, 21 min, 0.0 sec (37.350°) 

LON = 76° 0.0 min, 0.0 sec (76.000°) 

ALT = 20 ft 

TYPE = H, FLAG1 , FLAG2 

NORFOLK VOR/DME 

IDENT = ORE, MVAR = 7.5°W 

116.9 (chll6) frequency 

LAT = 36°, 53 min, 40 sec (36.894°) 

LON = 76°, 12 min, 0.0 sec (76.200°) 

ALT = 27 ft 

TYPE = T, FLAG1 , FLAG2 

FRANKLIN VOR/DME 

IDENT = FKN, MVAR = 6.75°W 

110.6 (ch43) frequency 

LAT = 36°, 42 min, 50 sec (36.714°) 

LON = 77°, 0.0 min, 30 sec (77.008°) 

ALT = 37 ft 

TYPE = L, FLAG1 , FLAG2 
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TABLE 3.- AIRCRAFT POSITIONS FOR TEST CASES 


A/C position #1 a,b 
LAT = 36°, 30s (36.05°) 

LON = 76°, 30s (76.05°) 

ALT = (a) 30,000 ft; (b) 15,000 ft 

A/C position #2 a,b 
LAT = 37°, 50s (37.0833) 

LON = 77° (77.0°) 

ALT = (a) 15,000 ft; (b) 10,000 ft 

A/C position #3 a,b 
LAT = 37° (37.0°) 

LON = 76° (76.0°) 

ALT = (a) 500 ft; (b) 10,000 ft 

A/C position #4 

LAT = 37° 31 min, 30 sec (37.525) 
LON = 77° 48 min, 0 sec (77.800) 
ALT = 30,000 ft 


f 
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INPUT 


VOR/DME MODULE 


Intermediate Output/Input 
requirement 


Power 

frequency; 

i-1,2 


•P; 


FI ight 

Dynamics 

Model 


Lat, Lon, ALT 

V 

g 




Power 

check 


Data base 
query for 
fi i-1,2 


fi 

Lat, Long, ALT 
*"Type, ident, flag 1 
flag 2 


DATA 

BASE 


validity check 
Geometry 
Calculations 
Bias 5 noise 
MVAR 


Auxiliary 
Input Requirement 


Output 



Bearing (variable) 

0m 

slant range (variable) 


Validity (discrete) 


Ident (text) 


Figure 1; VOR/DME SIMULATION MODEL 






FIGURE 2, STANDARD HIGH 
ALTITUDE SERVICE VOLUME 
(refer to FIGURE 6 for attitudes 


below 1000 feet (305 m) 

.100 nmi (185 km) 


60,000 ft. 
(18,288 m)* 


,130 nmi (241 km) 


^ - — — — ' 


14,500 ft 
(4.420 m)" 


FIGURE 3. STANDARD LOW 


40 nmi ALTITUDE SERVICE VOLUME 


(74 km) 


SURE 4 . STANDARD TERMINAL 

SERVICE VOLUME 

(refer to FIGURE Sfor altitudes 
below 1000 feet (305 m) 

26 nmi (46 km) 

/ V 12.000 ft 
(V J 1 (3.658 m) 


(refer to FIGURE 6 for altitudes 
below 1000 feet (305 m) 

40 nmi (74 km) 



NOTE: All tltvationi ihown art with ratpact 
to tha station's slta alavation (AGL). 
Matric Maasuramants ara givan for 
convanianca and ara approximations. 


— V . 


1,000 ft 






IN FEET 


FIGURE 5 . DEFINITION OF THE LOWER EDGE OF THE STA NDARD T 
(TERMINAL) SERVICE VOLUME 


H 


IBM 


■miiiiMiii 




!■■■■■■■! 


!■■■■■■■! 


■■■■■■■■■■I 


■■■■p^il 


!■■■■■■■■■■■■■■■■■■■■ 
■■■■■■■■■■■■■■■■■■all 
!■■■■■■■■■■■■■■■■■ v 


!■»■■■■■■■■■■■■■ ■^•■■■■I 

bas&ss:^ssss:::::s:i 


^^^■■■■■■■■■28555555555521 

■■■■■■■■■■■’’^■■■■■SSSSSSSS! 


■■■■■■■■■■■■858855P52SS5SSSSSSS5 




Immrnmm 

!!!S»!B!!S!SB!SS!!SSI: 


■■■■■■■■■■■■■■■■555! 


:Ss:s:s 25 s::s;::sis:ssass 5 S 5 sssH! 


DISTANCE TO THE STATION IN nml 

1 1 1 1 1 ” 

15 20 25 30 35 

DISTANCE TO THE STATION IN km 


FIGURE ft. DEFINITION OF THE LOWER EDGE OF THE STANDARD H 
(HIGH) AND L (LOW) SERVICE VOLUMES 


!■■■■■■■■■! 


■■■■■■■■■■■ 




■■■i 


IMI 


BISSK 

SB: n:s 

iaattsssassss::l 


!■■■■■■■■■■■■■■■■■■■ 




■■■■■^■■■■■SSSSSl 


!■■■■ 


■■■■■■■■■■■■■2285555 


iiiiiiffliKjjjiaijjjiitas:; 


as:::::: 


nmmmmuwumm 


!■■■■■■■■■ 


■■■■■■■■^■■■SS! 




■■■■■ 


•■■■■■■■■■■■■■85555! 


DISTANCE TO THE STATION IN nml 


DISTANCE TO THE STATION IN km 


IN METERS 






APPENDIX A 

A/C Power Failure to VOR #1 
Case, la (noise sources on) 

Input 

PVOR1 = false 
PV0R2 = true 
PDME1 = true 
PDME2 = true 
V0RF1 = 113.3 
VORF2 = 112.2 
A/C position = la 
Vy = 250 

Output 

Identl = blank 

Ident2 = CCV 

VFL01 = false 

VFL02 = true 

DFLOl = true 

DFL02 = true 

B01 = 0 +n v 

B02 = + n v + 213.3241 

DME01 = + N d + 88.0158 

DME02 = + N 0 + 56.6075 

Note 

n v and Nq terms should be examined to insure proper bounds and colored noise 
relationship. 
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APPENDIX A Continued 
A/C Power Failure to VOR # 2 
Case, lb (noise sources on) 

Input 

PV0R1 = true 
PV0R2 = false 
PDME1 = true 
PDME2 = true 
V0RF1 = 113.3 
VORF2 = 112.2 
A/C position = lb 
Vg = 250 

Output 

Identl = FAK 

Ident2 = blank 

VFL01 = true 

VFL02 = false 

DFL01 = true 

DFL02 = true 

B01 = + n y + 140.1943 

B02 = 0 + n v 

DME01 = + N d + 87.5665 

DME02 = + N q + 56.4257 

Note 

n v and N Q terms should be examined to insure proper bounds and colored 
noise relationship. 
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APPENDIX A Continued 
A/C Power Failure to DME #1 
Case, lc (noise sources on) 


Input 


PV0R1 = 

true 

PV0R2 = 

true 

PDME1 = 

fal se 

V0RF1 = 

112.2 

V0RF2 = 

110.6 

A/C posi 

tion 

Vg = 250 



Output 

Identl = CCV 

Ident2 = FKN 

VFL01 = true 

VFL02 = true 

DFLOl = false 

DFL02 = true 

1501 = + n v + 213.3241 

(502 = + n v + 124.3151 

0ME01 = 0 + N 0 

DME02 = + N d + 27.7270 

Note 


n y and N^ terms should be examined to insure proper bounds and colored 
noise relationship. 



APPENDIX A Continued 
A/C Power failure to DME #2 
Case, Id (noise sources on) 


Input 

PV0R1 = true 
PV0R2 = true 
PDME1 = true 
PDME2 = false 
V0RF1 = 116.9 
V0RF2 = 112.2 
A/C position = 3b 
Vg = 150 

Output 

Identl = ORD 

Ident2 = CCV 

VFL01 = true 

VFL02 = true 

DFL01 = true 

DFL02 = false 

B01 = + n v + 63.8901 

B02 = + n y + 187.7436 

DMEOl = + N d + 11.6823 

DME02 = 0 + N d 

Note 

Recheck noise terms for bounds and correctness of colored terms 
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APPENDIX A Continued 
Case 2, Station bearing checks 
(noise sources off) 

Input 

PV0R1 = true 
PV0R2 = true 
PDME1 = true 
PDME2 = true 
V0RF1 = 113.3 
V0RF2 = 112.2 
A/C position = 2a 
Vg = 200 


Output 


Identl 

= FAK 

Ident2 

= CCV 

VFL01 

= true 

VFL02 

= true 

DFL01 

= true 

DFL02 

= true 

B01 = 

70.9864 

B02 = 

309.6803 

DME01 

= 43.1044 


DME02 = 55.7343 



APPENDIX A Continued 
Case 3a, Valid geometry check 
(noise sources on) 

Input 

PV0R1 = true 
PV0R2 = true 
PDME1 = true 
PDME2 = true 
V0RF1 = 113.3 
V0RF2 = 112.2 
A/C position = 4 
Vg = 350 

Output 

Identl = Blank 

Ident2 = CCV 

VFL01 = false 

VFL02 = true 

DFL01 = false 

DFL02 = true 

B01 = 0 + n y 

B02 = + n v + 285.5277 

DME01 = 0 + N d 

DME02 = + N d + 86.5889 

Notes 

1) A/C inside cone of confusion for VOR/DME #1 

2) check noise terms 
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APPENDIX A Continued 
Case 3b, Valid geometry check 
(noise sources on) 


Input 

PV0R1 = true 
PV0R2 = true 
PDME1 = true 
PDME2 = true 
V0RF1 = 110.6 
V0RF2 = 113.3 
A/C position = 3a 
Vg = 120 

Output 

Identl = Blank 
Ident2 = Blank 
VFL01 = false 
VFL02 - false 
DFL01 = false 
DFL02 = false 
B01 = 0 + n v 
B02 = 0 + n v 
DME01 = 0 + N d 
DME02 = 0 + N d 



APPENDIX A Continued 
Case 3c, Valid geometry check 
(noise sources on) 


Input 

PV0R1 = true 
PV0R2 = true 
PDME1 = true 
PDME2 = true 
V0RF1 = 112.2 
V0RF2 = 116.9 
A/C position = 3a 
Vg = 120 

Output 

Identl = CCV 

Ident2 = ORF 

VFL01 = true 

VFL02 = true 

DFL01 = true 

DFL02 = true 

B01 = + n v + 187.7436 

802 = + n y + 63.8901 

DME01 = + N D + 21.0594 

DME02 = + N d + 11.5082 


Note 

check noise terms 
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APPENDIX A Continued 
Case 4a, Station failure check 
(noise source on) 


Input 

PV0R1 = true 
PV0R2 = true 
PDME1 = true 
PDME2 = true 
V0RF1 = 113.3 
V0RF2 = 112.2 
A/C position = 2a 
Vg = 175 

Output 

Identl = Blank 
Ident2 - CCV 
VFL01 = false 
VFL02 - true 
DFL01 = true 
DFL02 = true 
B01 = 0 

B02 = + n v + 309.6B03 
DME01 = + N q + 43.6810 
DME02 = + N d + 55.7343 

Note 


Data Base Flagll = false 



APPENDIX A Concluded 
Case 4b, Station failure check 
(noise sources off) 


Input 

PV0R1 = true 
PV0R2 = true 
PDME1 = true 
PDME2 = true 
V0RF1 = 113.3 
V0RF2 = 112.2 
A/C position = 2a 
Vg = 175 

Output 

Identl = FAK 
Ident2 = CCV 
VFL01 = true 
VFL02 = true 
DFL01 = true 
DFL02 = false 
B01 = 70.9864 
B02 = 309.6803 
DME01 = 43.8248 
DME02 = 0 

Note 

Data Base Flag22 = false 
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APPENDIX B 


Gauss-Marko Process for VOR/DME Error Models to Support ACS Project 

VOR: 

The VOR error model as synthesized previously is: 



Where: 


W v (t) = Gaussian White Noise with autocorrelation given as 


2(3 

R w (x) ^ e |w (t) w (t + %)] = vTT“ 6(x) 

V L J C G 


(Bl) 


O 

and f) is course roughness variance (specified) and 
W c = course roughness frequency (specified) 

Vq = aircraft ground speed 
N c (t) = course roughness error 

N r = a random constant defining receiver bias (a constant for each a/c but 
could change from run-to-run) 

i.e., N p = N(0, o p) with o p specified. 

N t = a random constant for each VOR station defining transmitter biases 

i.e., N t = N(0 , o t ) with o t specified. 

0 (t) = correct magnetic bearing from the VOR station to the aircraft. 

0 m = magnetic variation 

0 (t) = indicated VOR bearing 


C 


a(x )dx = 1. 
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APPENDIX B Continued 


The course roughness error, N c (t), is treated as colored noise 
correlation times inversely proportional to A/C ground speed 1 
autocorrelation of N c is given by 


R N (x) ' E[N c (t)N c (t +t)] = a 2 e 


9 W V 

= J- 0 c 9 


The power spectral density (W) and R w (x) are related bv 
transform pair, therefore, c c 


$ 


N ( w ) = | / r n (t ) cos wt dx 

C 71 J Q \ 


(t) = / $ (w) cos wr dw. 

y 0 N c 


With t = 0 equation (B4) reduces to 


r n (0) 


oo 

= f (w) dw = o ? c 


c 0 C 

Evaluating equation (B3) using equation (B2) gives: 

2 'VgM 


Y"> = l l \ e 


dx 


Therefore , 


o WJ r 9 

$ N (w) = I - ? ■ ? o a 2 

N c W 2 + W 2 V 2 c 

C b 


Thus N c can be simulated by computing the following response: 


n(t) 


J 2 “ cV °c 


V * S + « C V G 



-N c (t) 


where n(t) is Gaussian white noise with E(n) = 0; VAR(n) = 


wi th 
G* The 

(BZ) 

a Fourier 

(B3) 

(B4) 


(B5) 

(B6) 
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APPENDIX B Continued 


Digital Redesign of Continuous Process 

* A low pass digital filter that will solve eg. (B6) is given by: 


4 


N C (K) = «N c (k - 1) + (1 - a) K n(k - 1) 


(B7) 


-W r V r h , 

where a = e ^ (h = computation interval) and K is a constant gain to 
be determined which will insure the proper statistics for N c (k); and n(k) 
is a Gaussian white noise sequence with power spectral density 


0 


h/2u 



-n/h 


n /h 


i .e. , E(n(k) ) = 0 
VAR(n(k) ) = 1 


Therefore, 


^ (W) = h/2u u/h < W < u/h 


The Z-transform of (B7) is given as 


N C (Z) 
N(TT = 


( 1 . - «) K 

Z - a 


H(Z) 


(B8) 


Therefore, the power spectral density of the discrete sequence N c (k) is 


■» N (W) = | H(Z) | * n (W ) 


(B9) 


where by definition, Z = e 1 wh . The variance of N r is (using eqs. (B8) and 
(B9) 

u/h u/h 


L 


cj. . V W » dw ■ 

n /h c 


_ h C (1 - a 

T5T J .. . 2 

%/\\ 1 + a - 


- a) 2 K 2 dw 


2 a cos wh 


Integrating we obtain 


-2 _ 1 - a 2 

<J = i — T — K 
c 1 + a 
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APPENDIX B Continued 


We now determine K such that 
-2 2 

° c ~ a c (the desired variance based on definition of error model) 
Therefore , 



Substituting the above K into eq. ( 157 ) 

N c (k + 1) = , *N c (k) + /l - a ^ <? c n (k) 


The digital implementation of this difference 
random process for the course roughness error, 
bearing computation is given by 


equation simulates the proper 
Finally, the indicated VOR 


e (k + 1) = N (k + 1) + o(k + l) + n + n. + o 
c r t m 

where all quantities are defined above. 


DME ERROR MODEL 

The DME Error Model is synthesized as: 


W n (t) 



R D ( t ) 


Using the same techniques as 
the above provides the slant 
resulting equations used to simulate 


used in VOR development; the digital redesign of 
range computation as indicated by the DME. The 
this process are 




APPENDIX B Concluded 


where: 


N d (k + 1) =pN n (k) + /l - p 2 o R n(k) 
c c c 

P = e " h/x D; N b - N(0 >Od ) 


and 


E(n(k) ) = 0; VAR(n(k)) = 1 

. (Gaussian 

<T» (W) = 7T, -n/h < w < n/h \ + 

n 2,1 ( White 


The constants ap , tjp, and Tp are specified as before. 
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